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Abstract. Young, fast-rotating neutron stars are promising candidate sources for the pro¬ 
duction of ultrahigh energy cosmic rays (UHECRs). The interest in this model has recently 
been boosted by the latest chemical composition measurements of cosmic rays, that seem 
to show the presence of a heavy nuclear component at the highest energies. Neutrons stars, 
with their metal-rich surfaces, are potentially interesting sources of such nuclei, but some 
open issues remain: 1) is it possible to extract these nuclei from the star’s surface? 2) Do 
the nuclei survive the severe conditions present in the magnetosphere of the neutron star? 3) 
What happens to the surviving nuclei once they enter the wind that is launched outside the 
light cylinder? In this paper we address these issues in a quantitative way, proving that for 
the most reasonable range of neutron star surface temperatures (T < 10 7 K), a large fraction 
of heavy nuclei survive photo-disintegration losses. These processes, together with curvature 
losses and acceleration in the star’s electric potential, lead to injection of nuclei with a chem¬ 
ical composition that is mixed, even if only iron is extracted from the surface. We show that 
under certain conditions the chemical composition injected into the wind region is compatible 
with that required in previous work based on purely phenomenological arguments (typically 
~ 50% protons, ~ 30% CNO and ~ 20% Fe), and provides a reasonable explanation of the 
mass abundance inferred from ultra high energy data. 
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1 Introduction 

Magnetized and fast-spinning neutron stars were introduced early on as potential candidate 
sources of ultrahigh energy cosmic rays (UHECRs) [1—3]. The model has been revived recently 
[4-7], mainly because of new measurements of the chemical composition. These sources 
present many advantages compared to the more classical scenarios invoking gamma-ray bursts 
or Active Galactic Nuclei. First, the energy budget supplied by their population is relatively 
comfortable, due to their large rotational energy reservoir (E TOt 2 x 10 52 erg I Vi P_^, with I 
the star inertial momentum and P its spin period 1 , for an isolated new-born pulsar spinning 
close to the disruption limit), and their population density (h s ~ 3 x 10” 4 Mpc -3 yr _1 [8]). 
The energy injected into UHECRs is of order £uhecr ~ 0.5 x 10 45 erg Mpc 3 yr 1 [9], which 
implies that a fraction of order 10~ 4 of the neutron star population is required to achieve the 
UHECR flux level. In addition, the latest results reported by the Auger Observatory point 
towards a chemical composition of UHECRs that is not compatible with a light composition 
at the highest energies [10-12], The Telescope Array results seem to show the same trend 
within systematics [13-15]. Unlike AGNs and GRBs for which heavy nuclei production seems 
to be challenging [16-18], neutron stars, with their nretal-rich surfaces, offer a favorable site 
to produce heavy nuclei. 

One natural question is whether heavy nuclei stripped off the surface can survive the imme¬ 
diate environment of the star. This question was examined in Refs. [19-21] for the purpose 
of calculating the associated fluxes of neutrinos and gamma-rays, focusing on nuclei acceler¬ 
ated up to 10 15 - 16 eV in known local pulsars such as Crab and Vela. These authors showed 
that the thermal radiation field emanating from the neutron star surface was an issue for the 

1 Here and in what follows, quantities are labelled Q x = Q/ 10* in cgs units unless specified otherwise. 
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survival of high energy nuclei for surface temperatures T > 10' K. They also examined the 
effect of the non-thermal soft photon background on 10 15 ~ 16 eV nuclei, in the specific case of 
particle acceleration in the pulsar outer gap, at the light cylinder. 

Farther out, for distances much larger than the light cylinder, and for energies E > 10 18 eV, 
it was speculated that the curvature radiation losses might be avoided [1, 3, 6], and that the 
effects of the radiation field due to the nebula may not be such as to inhibit the escape of 
accelerated iron nuclei at the highest energies [6, 22, 23]. More important are usually the 
hadronic interactions in the supernova ejecta: Refs. [6] and [7] showed that even when these 
are dominant, escape of the highest energy particles is still allowed and a successful fit to the 
observed data can be obtained by tuning the initial composition (typically ~ 50% protons, 
~ 30% CNO and ~ 20% Fe). 

In this study, we re-examine the fate of heavy nuclei in the neutron-star immediate envi¬ 
ronment, before they are injected in the wind for acceleration to ultrahigh energies. One of 
our objectives is to estimate the chemical composition that is produced by photo-hadronic 
interactions and channelled into the wind, when pure iron is extracted from the surface of 
the star. 

To this purpose we calculate the photo-hadronic processes when iron nuclei are injected in 
the pulsar thermal radiation background. In agreement with Ref. [19], we find that for the 
most reasonable range of neutron star surface temperatures (T < 10' K), a large fraction 
of heavy nuclei survives the losses on the radiation environment of the star. Moreover, 
photo-disintegration leads to the production of secondary particles (intermediate nuclei and 
nucleons) that may help achieving a successful fit to the UHECR data, as required in Ref. [7] . 
The nuclei become then part of the pulsar wind, which is where most of the rotational energy 
lost by the star is thought to end up. The mechanism through which this wind is accelerated 
to its terminal Lorentz factor is not clear, but, from extrapolation of what we infer in young 
(though not newborn) pulsars, we assume that most of the wind energy is in the form 
of particle kinetic energy, rather than Poynting flux. For reasonable values of the pulsar 
multiplicity, the wind Lorentz factor is still high enough so as to guarantee that the heaviest 
elements have ultrahigh energy. We concentrate here on pulsars with rotation periods at birth 
P < 10 ms and mild dipolar magnetic fields (B ~ 10 12 ~ 13 G), as this parameter space was 
found to be ideal for successful acceleration of UHECRs and their escape from the surrounding 
supernova ejecta [6]. The existence of such pulsars is supported by pulsar population studies 
[24], and could be confirmed by the observation of X-rays and gamma-ray signatures from 
ultra-luminous supernovae [25-27]. Some authors also suggest that the radio-emission from 
the Crab nebula imply that the Crab pulsar was born with a 5 ms period [28]. 

The paper is organized as follows: in Section 2 we discuss the processes through which 
particles can be extracted and accelerated within the light cylinder of a pulsar. Taking into 
account the various radiation losses at play, we calculate the maximum Lorentz factor that 
particles can reach, mainly determined by the balance between acceleration in the electric 
field of the pulsar and curvature losses (this Lorentz factor is important when compared with 
the threshold of photo-hadronic interactions). In Section 3, we detail our calculations of 
the production of secondary products of photo-disintegration, that determine the chemical 
composition of particle injected in the wind. Section 4 sketches the later acceleration of 
particles in the wind, and presents the spectra that will be injected in the nebula and/or 
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the supernova ejecta region for further processing. We discuss our results and conclude in 
Section 5. 


2 Ion extraction and acceleration within the light cylinder 

Neutron-star surfaces are believed to be composed of anisotropic, tightly bound condensed 
matter. Calculations show that the outer layers should be composed mainly of long molecular 
chains with axes parallel to the magnetic field [29, 30] . The chains are thought to be composed 
of 56 Fe ions forming a one-dimensional lattice with an outer sheath of electrons. The binding 
energy of the ions can be estimated as ~ 14keV and the lattice spacing as l rs_/ 10 9 cm [31]. 
Ions can thus be stripped off the surface provided that they experience a surface electric held 
of £ 0 = 14keV/ (Zel) ~ 5.4 x 10 11 Z^lz\ Vcm" 1 . 

However, ions can also be stripped off with milder electric fields if the surface is bombarded 
by particles (for instance leptons accelerated towards the star surface after the cascading 
processes in the magnetosphere are started), and/or if they can be boiled off the surface by 
stellar heat [31, 32], In the case of a pulsar with millisecond periods at birth, the electric 
held that can be provided at the surface of the star can be estimated as (see, e.g., [3]) 

£ = 27T p* B ~ 6.3 x 10 14 B^R+fiPzl V cm -1 , (2.1) 

where B is the dipole magnetic held strength of the star, P* its radius and P its rotation 
period. In practice, this electric held does not persist on large distances, as it is expected to 
be readily screened out by pair cascading. 

However the total potential drop that is available, at least in principle, for particle accelera¬ 
tion in the magnetosphere is 2 : 

$ = 2n p2 c ? ~ 6 ' 6 >< 1019 B 13 R l,6 P -3 V , (2.2) 

which corresponds to a maximum achievable particle Lorentz factor: 

^=^*= 7xi ° 10 ! b - p -t (2 - 3 > 

In reality, the maximum energy achievable by ions within the corotating region will be 
limited by losses, the most severe among which is curvature loss. Assuming that the to¬ 
tal potential drop <f> is available for particle acceleration over a gap of length with 
= cP/ (27r) ~ 4.8 x 10 6 P__3 cm the light cylinder radius, the equation describing the 
particle energy evolution, subject to acceleration and curvature losses can be written as: 

±L = _Ze*_ 2* _ _Sn 

dt ArripC 2 £ P 3 cP 2 Am p c 2 

throughout this work, we have adopted the convention used in Ref. [3] for the geometric pre-factors of 
the electromagnetic energy loss rate. See in particular its footnote 5. 
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The maximum Lorentz factor achievable for radiation limited acceleration is then found 
equating the two terms in the r.h.s. of Eq. 2.4, and turns out to be: 


7curv = (H||) ' ~ 1-1 X 10 8 Z 2 - 1/4 r 1/4 < 4 ^-3 1/4 <6 • (2-5) 

The actual maximum energy that particles can reach at any time within the corotating 
magnetosphere will be set by: 

7max = min( 7 fpd , 7 curv ) . ( 2 . 6 ) 


One thing to notice in this expression is that 7 CU rv depends very weakly on the fraction of R l 
over which the gap extends (£). On the other hand, the acceleration time to a given energy, 
t acc ( 7 ), has a much stronger dependence on the unknown 0 


t&cci'j) 
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= 5 x 10~ 9 s 


7 


A, 


56 


10 6 Z26 


B 


-1 

13 


Kl 


p-3 a 


(2.7) 


On the other hand, while in general £ is rather uncertain, depending on pair creation in the 
polar cap, the constraint < £ < 1 is particularly tight for the fast-spinning pulsars we 

are concerned with: in this case R+/Rl ~ 0.2 R^qP^ 3 implies £ = 0 ( 1 ) and that the gap 
cannot be far from the surface of the star. 


In this paper, we will work under the most classical assumption that particles are accelerated 
close to the stellar surface. For scenarios in which the acceleration happens at the light 
cylinder or further out, the thermal background could be too anisotropic and diluted to play 
a role. On the other hand, for the outer gap scenario [33, 34], a non-thermal soft photon 
background could exist that could play a similar target role as the thermal background 
[20, 21]. Due to the even larger amount of uncertainties it entails, we will not consider this 
scenario further in the paper. 

One of the key points we want to make in this paper is to stress the crucial role of the 
Lorentz factor of the wind: while in the literature one can often find the argument that 
pulsars cannot be sources of CRs with Lorentz factor above 7 max (Eq. 2.6), because of the 
dramatic curvature losses inside the light cylinder, we argue that if 7 max < 7w (where y w is 
the Lorentz factor of the wind, that we express later in Eq. 4.1), once particles end up in 
the wind, they get advected with it at the Lorentz factor y w , irrespective of the energy they 
reached in the magnetosphere. We will discuss the implications of this point in Section 4. 

Notice also that while photo-disintegration of nuclei may become very important inside the 
light cylinder, these interactions conserve the Lorentz factor of the nuclei involved. Photo- 
pion production of protons can in principle change the Lorentz factor of individual nucleons, 
but this process is very slow and its effect is found to be negligible. 


3 Photo-hadronic interactions on the thermal radiation of the star 

Ultrahigh energy ions can experience photo-pion production or photo-disintegration in the 
thermal radiation fields generated by the star. As we will see in the following, the latter 
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process is much more important than the former and leads to the production of progressively 
lighter nuclei, with interactions that change the number of nucleons while conserving the 
particle Lorentz factor. Before discussing in detail the consequences of this process on the 
composition of the UHECR flux from a newborn pulsar, we first recall briefly our current 
knowledge on the temperature of the pulsar surface at early stages, and estimate the key 
quantities involved in the photo-hadronic processes. 


3.1 Photo-hadronic interaction timescales 

The thermal structure of the neutron star upper layers is not well understood, and constraints 
from observations are scarce (e.g., [35-37]). Neutron stars with magnetized envelopes made 
of accreted material have local surface temperatures of typically T 1q5.5 i anc [ £ or 
magnetized objects, show a profile that can vary as a function of the polar angle. 

The thermal photon energy density, 17th, can be expressed as a function of the neutron star 
surface temperature T as 17th = aT 4 , a being the radiation constant. This thermal component 
peaks at energy e 7;t h = ksT ~ 86 eV T§. 

One can model the photo-disintegration process for nuclei with a delta function (representing 
the Giant Dipole Resonance with a threshold energy of order Za^ = 18.3l7l5 6 0 ' 21 MeV and 
a width of Ae = 8MeV), and a roughly constant tail, with a factor of 8 times lower cross 
section, at higher energies [38, 39]. The cross-section at the resonance can be written as o a~< ~ 
8 x 1(T 26 R 56 cm 2 . The threshold energy is of order (ei + 62)/2 = e^ 7 = 18.3lA^ 6 °' 21 MeV and 
62 — ei = Ae [38-40]. The corresponding threshold Lorentz factor for particles thus reads 

7Tthre S = —~2.1x 10 5 R 5 - 6 0 - 21 Tf 1 . (3.1) 

e 7,tli 

For typical parameters appropriate for a highly magnetised fast spinning pulsar, this thresh¬ 
old energy is lower than the maximum energy 7 max reached by particles accelerated in the 
gap (Eq. 2.6), implying that most particles are in a condition to experience photo-hadronic 
processes. 


The photo-hadronic interaction time on the neutron star thermal background can be ex¬ 
pressed for a nucleus of mass number A, with Lorentz factor 7 a and velocity j3a, as 


tA 


c [°° de dn(e,T) f 1 a A -y(e')(l-Mdn 

- ^min ^ e 


(3.2) 


with e 1 = 7a( 1 — /?A/^)e (primed quantities are in the nucleus rest frame). Here fj, = cos 6, 
with 0 the angle between the direction of the cosmic ray and the photon (in the lab frame). 
The maximum angle from which a thermal photon produced at the surface of the star can 
hit a nucleus at a distance ct from the surface is then (t) = [1 — R 2 /(R* + ct) 2 ] 1 / 2 . The 
anisotropic orientation of the photon background field radiated from the stellar surface is 
taken into account by a suitable choice of the integration limits in the integral over //. The 
photon number density dn(e, T)/de ~ L r th/ e 7 ,th ( l( e — e 7 ,th) is calculated at the surface of the 
star. 
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Since we will mostly be concerned with non resonant interactions, the relevant interaction 
timescale can be approximated by 


tA = t* A XA{t) , 


(3.3) 


where 


t*A = 


O'A'y Ufa 


^ 7 ,th 


-1 


6 x 10 0 s A 5 g 1 T 6 3 


and 


(t) -! = I (! “ Pah) dh if ei < 7^(1 - Pah)* < e 2 
l £«(() C 1 ” d h- if 7^(f - Pah)(- > £2 • 


(3.4) 


(3.5) 


We define the instantaneous opacity of the interaction as t a = t/t A . For the sake of illustra¬ 
tion, at the stellar surface, where the radiation field can be considered as roughly isotropic 
(XA(t = 0) 1 for interactions out of the Giant Dipole resonance, second line in Eq. 3.5), 

this quantity reads r A = Ri./(ct A ) ~ 0.5 A^q Tg . 

For photo-pion production by nucleons, the same scheme as for the Giant Dipole resonance 
can be adopted with (ei + € 2 )/2 = 0.3 GeV, £2 — ei = 0.2 GeV, and a vl = 4 x 10~ 28 cm 2 [41]. 
The much smaller cross-section compared to photo disintegration makes this process irrele¬ 
vant in all cases, as we will see in the following section. 


3.2 Primary and secondary nuclei ratios 


In what follows, we calculate the number of primary and secondary nuclei that become part 
of the pulsar wind, after particles have been extracted from the surface of the star, and have 
suffered photo-disintegration. We assume for simplicity that nuclei lose one nucleon in each 
interaction. We assume that pure iron is injected at time t = 0 at the surface of the star. 


1 

1- 

■HO 

l~tk 

lt =0 


The evolution of the iron abundance then follows: 

Afe = N Fe (t = 0) exp 


For the other species with mass number A > 1, one can write 

d-Wi Na _ N a+ i 

df t A t A+ i 


which, assuming N A (t = 0) = 0 for A < 56, leads to: 
N a = exp 


1 /■* dt' 1 

f l n a+ 1 

’ 1 /■*' dt" ’ 

[ t* A Jo x(t')\ 

/ £ 6XP 
Jo T4+1 

t* A Jo X(t")_ 


dt' . 


(3.6) 


(3.7) 


(3.8) 


Figure 1 shows the evolution of N A calculated via this equation with t A + and x defined as 
in Eqs. 3.4 and 3.5. For surface temperatures exceeding 10 6 K, photo-disintegration starts to 
be effective and with increasing temperature an increasingly large fraction of the initial iron 
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Figure 1 . Composition of UHECR nuclei after photo-disintegration in the thermal radiation field of 
the star (solid lines). Extraction of pure iron at t = 0 is assumed. Eqs. 3.6 and 3.8 are solved for 7 (f) 
determined by Eq. 2.4 (with £ = 1) up to the value 7 ma x set by Eq. 2.6. The fractions of different 
species are represented by solid lines, with colours as specified in the legenda of each panel, while 
dashed lines refer to the opacities ta = t/tA • From left to right, top to bottom, temperatures are 
T = [1,2,5,10] x 10® K. The vertical line indicates the time at which particles reach the light cylinder 
i?LC -1 and the dotted horizontal line indicates Na = 1 . 


is transformed into lighter nuclei. When the star temperature is of order 10'K, no iron is 
left at the light cylinder, nor elements with Z > 10. In this case the composition is mostly 
protons and ~ 10% CNO. The trend shown by the opacity of the different nuclei is readily 
understood: initially it simply increases with time, but then it saturates and starts dropping 
due to the decreasing photon held density. The maximum opacity is reached at a time which 
is independent of nuclear species and stellar temperature, as can be readily understood from 
Eqs. 3.4—3.5. 

One thing to notice is that, as we already anticipated at the end of the previous section, the 
photo-pion process is always irrelevant, as shown by the dashed black line representing the 
corresponding opacity in the different panels. 

These photo-hadronic interaction calculations are done under rough approximations. We 
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have considered a cross-section composed of a delta function (to mimic the Giant Dipole 
Resonance) and a plateau at high energies (for higher-energy processes), and assumed that 
only one nucleon is liberated in each interaction. In practice, at the energies that we consider, 
we barely fall in the Giant Dipole Resonance regime, but in the “plateau”, where Quasi- 
Deuteron and Baryonic Resonance processes dominate. The products stripped from the 
nuclei via these interactions are likely not single nucleons, but deuterons or alpha particles. 
Fragmentation of the nuclei in several lighter nuclei is also possible. This constitutes a limit 
of our calculation. A more thorough numerical calculation with tabulated processes would 
enable to better assess the resulting composition. A better treatment, however, would not 
change our main result: heavy nuclei would still survive for temperatures T < 10' K. 


4 Spectra injected in the wind 


Particles accelerated in the gap are injected in the pulsar wind, and reach higher energies, 
with fluxes that we estimate below. They subsequently reach the nebula region and the 
supernova ejecta, where they encounter other radiative and baryonic fields. As mentioned in 
the introduction, the effects of these backgrounds were calculated extensively in [6, 7, 22]. 
More precisely, Ref. [7] estimated that the typical composition that has to be accelerated 
in the wind region before being altered by the supernova ejecta and the propagation in the 
intergalactic medium was ~ 50% protons, ~ 30% CNO and ~ 20% Fe. These authors discuss 
that this result is rather robust (the composition weakly depends on the supernova parameters 
or on the narrow range of pulsar parameters allowed for successful UHECR acceleration). 

Neutron stars spin down and their rotational energy is channelled via their winds towards the 
outer medium. In the following we assume, as the observations suggest for the Crab pulsar, 
that the electromagnetic luminosity, E p = (8n B 2 R^) / (3c 3 P 4 ) [42], is efficiently converted into 
kinetic luminosity, N^mc 2 . We write the latter as: N^mc 2 = 2Km e c 2 A r GJ7w (1 + xa), 
with IVgj = e~ 1 (E p c) 1 / 2 the Goldreich-Julian charge density [43], and xa = btla/ (2nZm e ), 
where tha = Am p is the ion mass and k the pair multiplicity, namely the number of pairs 
produced for each electrons that leaves the star. 

We can then write the wind Lorentz factor as [22]: 


7w (t) 




(1 - Vb) 


4 

Nmc 2 


~ 3 X 10 9 (1 - r] B )(l + X A ) ‘(1 + t/t sd ) 1 ^ P i,- 3 B ^ R l ,6 


(4.1) 


Here, we have expressed the pulsar spin down luminosity E p (t ) = E p o/(l + t/t s d) 2 (as¬ 
suming a breaking index of 3), where the initial luminosity is E p q = E TO t/t sd — 6.4 x 
10 44 P~\B 2 13 R^ 6 erg/s, and the spin-down time is 


^sd 


9 Ic 3 P 2 
8vr 2 B 2 R e 


~ 3.1 x 1() 7 s R P By 2 R,^ gT ) 2 _3 . 


(4.2) 


t]b represents the fraction of magnetic energy injected into the energy conversion region. 
Observations indicate that close to the pulsar wind nebula t/b <C 1 [44]. 
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Depending on the value of k, that can range between 10 — 10 8 in theory (a highly debated 
quantity, see e.g., [44]), the energy conversion can be efficient enough to enable iron to reach 
energies at neutron-star birth (t = 0 ) 

E 0 ~ 2.3 x 10 20 eVA 56 r] P^ 3 Bi 3 Rl fi , (4.3) 

assuming xa < 1 (we have, e.g., XFe ~ 0.2 for k = 10 4 ) and ijb <C 1. To take into account all 
the uncertainties on the luminosity conversion efficiency, we have introduced the efficiency 
factor Tj < 1 . 


As the neutron-star spins down, cosmic rays flowing with the wind will have an energy 

EcR,{t) = E 0 (1+t/tsd) 1 (4.4) 

~ 1.2 x 10 20 eV R~q £ 7 . 5 “* for t > t s a- 

Channelling the Goldreich-Julian charge density into particles and taking into account the 
neutron-star spin down rate (neglecting gravitational wave losses), one can write the cosrnic- 
ray injection flux as: 

= J™ dtN GJ (t)6 (E - E CR (t )) = Ncj( £ )tsd . (4.5) 

In order to compute the spectrum and composition emerging from the star, we can use the 
results of the previous section and Eq. 4.5 as long as a number of conditions are satisfied: 


• the star surface temperature should not decrease appreciably below the value initially 
assumed 

• the electric field at the star surface, which decreases with time as P ” 1 must stay larger 
than the value required for extraction of iron nuclei (which is also changing and likely 
increasing with time, due to the decreasing temperature) 

• the acceleration time (Eq. 2.7) needed to reach qth, which has a strong dependence on 
the pulsar period, oc P 3 , should still be short enough so that 7th is reached before the 
time corresponding to the peak in the opacity (which always occurs at t « 10~ 2 Pl(0)/c, 
as one can see from Fig. 1). 


While the evolution of the star surface temperature is not very well constrained, it seems 
reasonable to assume that this stays constant for the first few xlOOyr. If this is the case, 
then all other conditions are satisfied for t < 100 yr, and we decided to limit ourselves to this 
time. This causes us to be able to reliably compute the spectrum of iron only down to an 
energy E ~ 10 18 eV (see Eq. 4.5). This does not mean that the pulsar model predicts a low 
energy cut-off for the spectrum of iron at around this energy, with the UHECR composition 
predicted by this model becoming progressively lighter at lower energy. A low energy cut-off 
for iron is expected, due to the fact that condition 2 will be violated at some point, but 
currently we are simply not able to assess this matter, due to the large number of unknowns. 
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Figure 2. Cosmic-ray spectra EdNA/dE and composition estimated at time t w after injection with 
t w = 100 x Rl/c and an injection time < 100 yr. This limit on £; n j, which reflects in the iron spectrum 
to be cut off at ss 10 18 eV, was chosen because at later times the unknowns are such as to make our 
calculation unreliable (see text for details). The results refer to a case in which only iron is extracted 
from the star surface and photo disintegration is taken into account as described in Section 3. From 
left to right, top to bottom, temperatures are T = [1,2, 5,10] x 10 6 K. 


Figure 2 shows the spectra of the primary and secondary nuclei produced in the wind at 
times f w = 100 x R\,jc after injection, for a primary iron nucleus injection at t = fi n j with 
t - m j < 100 yr. The time t w is chosen to be much longer than the photodisintegration loss 
time scale, to ensure that these spectra reflect the composition that will be injected in the 
wind and leave the star with the wind Lorentz factor. In our calculations we use the output 
obtained in Section 3 and Eq. (4.5) normalized to the total number of injected nuclei: 


cEVa _ N a f 7 

d E N Fe (t = 0) 6XP V 


(4.6) 


where y w is defined in Eq. (4.1) with xa < 1 for all nuclei. 


Figure 3 gives the fraction of each nuclear species as a function of the neutron star surface 
temperature T, assuming again that only iron is directly extracted from the stellar surface. 
This figure demonstrates that a significant fraction > 10% of heavy nuclei survive the photo- 
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disintegration process on the neutron star thermal background, for temperatures T < 5 x 
10 6 K. For higher temperatures, lighter nuclei are injected together with protons at a rate of 
~ 5%. It is interesting to note that this is a reasonable range of temperatures based on our 
knowledge of neutron star surface, and it would also provide a composition compatible with 
that required by Ref. [7] to fit the UHECR data. 



I [K] 

Figure 3. Fraction of each nuclear species injected in the wind at time t w for energies 7 < 7 , 4 . max, 
as a function of the neutron star surface temperature T, assuming initial pure iron injection. 


5 Discussion and conclusions 


The sources of UHECRs are still to be found. The observational situation has changed quite 
a bit in the last few years, leading to new challenges. The main new pieces of information 
that have arisen from Auger data can be summarised as follows: 1) The mass composition, 
as measured from the the elongation rate A max (E) and its fluctuations, drifts from a lighter 
one at E ~ 10 18 eV to a heavier one at the highest energies [7, 10, 15]; 2) The injection 
spectrum required to fit the fluxes and the mass composition at the same time seems to be 
harder than the typical ~ E~ 2 inferred for most astrophysical sources [45]; 3) If taken to 
scale as the charge Z of the nucleus, the maximum energy is bound to be ~ 5 x 10 18 Z eV 
[45], thereby relaxing the tight constraints on the maximum energy in case of pure proton 
composition. 

All of these three requirements appear to be fulfilled if fast spinning very young pulsars are the 
sources of UHECRs: as recognized already in early work on this scenario [2, 3], the spectrum 
of particles accelerated in the pulsar magnetosphere is very hard, typically oc E~ l . Since 
the particles that get energized are extracted from the the iron-rich surface of the star, it is 
likely that the initial composition of the accelerated particles is dominated by Fe nuclei [2], 
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although propagation in the ejecta of the parent supernova may change the abundances of 
lighter elements [7]. 

Within the pulsar scenario for the origin of UHECRs, one fundamental question that has 
never been convincingly addressed, in the relevant parameter range, is that of Fe survival in 
the hostile environment provided by the close magnetosphere. This has represented the main 
goal of the present paper: we have calculated the effects of the immediate pulsar environment 
on heavy nuclei, focusing on the dominant background, represented by the thermal radiation 
field from the star. The electric field due to the fast rotation of the neutron star gives rise 
to the extraction of iron nuclei from the star surface. Even in the simple case of pure Fe 
composition of the stripped off material, we found that a mixed composition is likely to be 
finally produced, as a result of photo-disintegration processes in the magnetosphere of the 
pulsar. For reasonable assumptions on the star temperature (T few x 10 6 K) and stellar 
parameters, it is possible to have a relative abundance of the injected elements that is in 
qualitative agreement with that required to reproduce the mixed composition observed by 
Auger (see Ref. [7]). In particular, the intermediate mass composition suggested by the latest 
Auger data [11, 12], can be easily accommodated for stellar surface temperatures around 
~ 5 x 10 6 K. 

Our calculations also show that at times of order ~ 100 years, it becomes challenging for the 
induced electric held at the star surface to overcome the binding energy of iron nuclei in the 
lattice, so that a low energy cutoff (or suppression) is induced in the spectra of injected iron, 
typically in the EeV energy range. Lower energy particles can still be produced as a result 
of photo disintegration of iron nuclei in the magnetosphere. 

Once the nuclei have been stripped off the surface of the star, they get accelerated in the 
gap, assumed to be distributed on a region of size Rl (radius of the light cylinder) from the 
stellar surface. While being accelerated, particles suffer photo-disintegration processes on 
the photon background emitted by the hot surface of the star, and lose energy to curvature 
radiation. Curvature energy would limit the Lorentz factor of the nuclei to ~ 10 8 , which for 
fast spinning neutron stars is appreciably lower than the value reachable by using the full 
potential drop. However, the nuclei that reach the light cylinder region eventually end up 
in the wind of electron-positrons propagating outwards. If the multiplicity of electrons and 
positrons is < 2 m p /m e , the maximum energy of iron nuclei (and of the secondary nuclei that 
are produced by photo disintegration) is in the range of interest for UHECRs. Notice that 
by construction the maximum energy of the secondary nuclei scale proportional to the mass 
number A instead of the charge Z of the nucleus. 
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